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ABSTRACT 


The report outlines progress of the research conducted during the period 
of June 1984-January 1985. During this period Items 2.1, 2.2 and 2.5 In the 
project proposal have 'jeen mostly accomplished. Moreover, a closure of 
third-order was pursued to Improve the diffusion rate of the Reynolds-stress 
transport equations for the prediction of recirculating and reattaching flows. 


r 


TABLE OF CONTENTS 

ABSTRACT 1 

TABLE OF CONTENTS 11 

NOMENCLATURE Ill 

1 . INTRODUCTION 1 

2. MATHEMATICAL FORMULATION 3 

2.1 Reynolds-Stress Equations 4 

2.2 Third-Order Closure Model 8 

2.3 Other Models of Third-Order Closure 10 

3. NUMERICAL MODEL 12 

3.1 Numerical Method , ... 12 

3.2 Boundary Conditions 13 

3.3 Computing Details 14 

4. RESULTS AND DISCUSSION 15 

4.1 Mean Velocity Profiles 15 

4.2 Reynolds Stresses 17 

4.3 Triple Velocity Products 18 

5. SUMMARIZING REMARKS 20 

6. WORK IN PROGRESS 22 

6.1 QUICK Method 22 

6.2 Consolidation of the RSM 24 

7. APPENDIX ..... 27 

7.1 Near-Wall Values for Reynolds Stresses 27 

8. REFERENCES 29 

9. FIGURES 31 


- 11 


t 


1 


B l» b 2* b 3 

C] , C2, C3, Cq, 
Cs* C c -j , C c 2, C, 

D u 

G 

H 

k 


P 

P 



Re o 


u 


U 



u 

T 


v 


V 


x 


x 


r 


y 

Y 


0 


z 


4 u 

AA 


c 


NOMENCLATURE 

constants used In turbulence model 
coefficients used In turbulence model 

diffusion of iTjTJJ 

generation rate of turbulence kinetic energy, k 
step height 

turbulence kinetic energy 
pressure fluctuation 
mean pressure 

generation rate of Reynolds stresses 
generation rate used in pressure-strain correlation 
Reynolds number based on momentum thickness, 0 
fluctuating velocity 
mean velocity 

mean velocity at the inlet 
friction velocity (=ST~7p) 

W 

fluctuating velocity in y-direction 

mean velocity in y-direction 

Cartesian coordinate 

length of reattachment 

Cartesian coordinate 

height of inlet flow section 

coordinate normal to the wall 

Kronecker delta 

area of numerical cell 

energy dissipation rate 
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Subscripts 
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grid expansion factor In x-direction 
dynamic molecular viscosity 
density 

Prandtl numbers for k and c 
wall shear stress 
pressure-strain correlation 


tensor notations 
wall values 
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1 . INTRODUCTION 


Flow In a channel which has steps or sharp bends occurs In many engineer- 
ing applications, such as airfoils with separation bubbles, combustors, heat 
exchangers, etc. While there have been significant contributions In this 
field, current understanding of the process Is still relatively poor, partly 
because the flow In these geometries is complex, and partly because theoreti- 
cal models developed to date are limited in predicting a wide range of para- 
meters In separating, reattaching, and recirculating flows. 

1 2 

Bradshaw and Wong and Smyth measured turbulence quantities in the 
reattaching and redeveloping regions and reported that the center part of the 
reattaching shear-layer still has most of the characteristics of a free shear 
layer as many as 50 step heights downstream of the reattachment point. This 
suggests that the large-scale eddies that develop in the separated free shear 
layer persist in the reattaching and redeveloping regions. However, the 
turbulent shear stress in the reattaching layer is much larger than the shear 
stress in an ordinary plane-mixing layer and it decreases rapidly near the 
reattachment point. 

3 

Eaton and Johnston focused measurement on the separated flow over the 
backward-facing step. A comparison with a plane-mixing layer showed that the 
normal components of Reynolds stresses in the shear layer near the reattach- 
ment point were 10-20% higher than typical values measured in plane-mixing 
layers. However, the values of shear stress were of about the same order as 
those of plane-mixing layers in most of the region upstream of the reattach- 
ment. It was also shown that the reattaching layer was similar to a plane 
mixing layer downstream of the reattachment region. These observations Indi- 
cate that many parameters could be predicted in the reattaching shear layer by 


using simple mixing-layer data; and, consequently, the turbulence models used 
In a mixing-layer flow would still be adequate for the reattaching shear layer. 

Theoretical studies for turbulence modeling for a wide range of flow 
patterns are available, Among the existing models, the so-called Reynolds- 
stress model of turbulence, which consists of transport equations for all the 
Reynolds stresses, provides Individual stress behaviors, while the viscosity 
based-models like the k ~ c model cannot account for such behaviors 
because an Isotropic assumption Is used. A model of the Reynolds stresses was 

4 

first proposed by Rotta , and has been developed and improved by a number of 
5-8 

researchers. The correlation of pressure-strain was proposed by Naot et 

C Q 

al . and Launder et al. Naot et al. evaluated the pressure-strain 

correlation term by Integrating over space after Inserting a conjectured form 

for the two-point correlation functions, whereas Launder et al, obtained the 

results by assuming a fourth-order tensor consisting uf linear Reynolds-stress 

elements. As a simpllfieu model of the Reynolds-stress closure model, an 

algebraic stress model which does not possess both convection and diffusion 

g 

terms was developed in a similar manner by Rodi. 

In the authors' previous paper^® several proposed closures for the 
Reynolds-stress model were tested for heat transfer characteristics along 
walls of axisymmetric sudden expansion pipes. After computation using several 
different models, it was found that the pressure-strain correlation proposed 

g 

by Launder et al. showed slightly better agreement with experimental data 
for Nusselt number distributions along the pipe walls. Moreover, an incor- 
poration of the wall correction terms In the pressure-strain also improved 
predictions by about 5 to 10 percent. 
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Although the Reynolds-stress closure model (RSM)* was formulated In the 
previous work 10 , the behavior of the stress components was not examined due 
to a lack of suitable experimental data on the Reynolds stresses. In this 
study the previous model was reviewed again, and Reynolds-stress levels 
computed by the model presently developed were compared with experimental data 
available in the literature. The results of the Reynolds-stress levels were 
also compared with the computations obtained by using the standard k ~ c 
model and the algebraic stress model. 

Meanwhile, the authors also noticed the fact that the triple velocity 
products In the separated shear layer begin to show a marked drop in the 
near-reattachment region due to a suppression of the large-scale eddies as the 
flow approaches a solid wall. This indicates that the third-order tensor 
quantities, which have direct Influence on the diffusive action of the 
Reynolds stresses, must be evaluated carefully for a reattaching shear layer. 
This study was added to the proposed research before the Reynolds stress model 
is extended further for heat transfer problems. In this way the diffusion 
rate of the second-order closure model of the Reynolds stresses will 
eventually be treated by the third-order closure model. The mathematical 
formulation of all the models is described in the following section. 

2. MATHEMATICAL FORMULATION 

The steady, two-dimensional form of continuity and momentum equations 
describing the flow field are used in this study. These are given as: 


*This is not the complete Reynolds-stress model because the stresses in the 
momentum equations were evaluated by using the Boussinesq viscosity model. 
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Continuity Equation: 
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Momentum Equation: 
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2.1 Revnolds-Stress Equations 

If the Reynolds st» esses, u^uj , were computed through the solution of 
the Reynolds stress transport equations, these values could be simply inserted 
Into Eq. (2); and the final flow field could be obtained by solving Eqs. (1) 
and (2) and the Reynolds-stress equations iteratively. However, considerable 
difficulty was encountered In solving the momentum equations due to the 
numerically unstable nature of these equations. Since the Reynolds stresses, 
puTuT, In Eq. (2) are evaluated through the source terms of the discre- 
tization equations, the diffusion term of the momentum equations contains only 
molecular diffusion rates. Consequently the convection rate of the momentum 
predominates over the diffusion rate, causing inherent Instability in the 
solution. This problem is overcome by bringing the Reynolds stresses into the 
diffusion term with a suitable modification. However, since this modification 
has been successfully carried out just recently and is still in the testing 
process for universal usage, the conventional Boussinesq effective viscosity 
was tentatively adopted to demonstrate the models developed so far. The 
description of the new diffusion part is given in the "WORK IN PROGRESS" 
section at the end of this report (Sec. 6.2). 
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The closure of the Reynolds stresses In £q. ( 2 ) was achieved as follows: 


0 i2 ®Uj 2 
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where k and c are obtained by solving the standard k ~ e model as shown 
below. 
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Here the generation rate of the turbulence kinetic energy can be expressed as 

au. 


6 ‘ • u i u j s? 


(6) 


The values for the constants used above are given as 


c 

M 

C cl 

C c2 

*k 

a 

c 

0.09 

1.44 

1.92 

1.00 

1.30 


The transport equations for the Reynolds stresses are given as 
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Equation (9) was approximated by the form given by Rotta 4 as: 

c 1j “ 3 6 1j c (12) 

4 

The pressure strain term (10) was determined combining Rovta's linear 
return-to-isotropy hypothesis and the linear approximation of Launder et 

g 

al. The results are given as: 
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Cj «s 1.5, C 2 » 0.4 
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and 



au k aU. 

-<“f u K 5^ + u j u k S^> 


(15) 


pressure-strain correlation with a near-wall correction* Is given as 

1 

♦u.w ■ t°- 125 k ( vj - f k *u> + 0015 ( p u - v 1 5? < 16 > 

where 2 Is a coordinate taken normal to the wall. Thus, If the point of 
Interest Is close to two plates (say z 1 to one plate and z ? to the other), 
then z Is chosen to be 

z « l/(l/z 1 + l/z 2 ) (17) 


The diffusion rate (11) is evaluated as: 


’u 


3 Vi 


“ C D ax,/ c u k u m ax 
k m 


(18) 


where Cp is a constant with a value of 0.25 which was determined by com- 
puter optimization. 

It should be noted that both convection and diffusion rates of the 
Reynolds stresses are neglected in the so-called algebraic stress model (ASM) 
whereas they are kept in the present RSM. 


^Although this correlation was used for the computation of the Reynolds 
stresses, the contribution of this term was relatively minor. That is why 
this term was dropped when the triple velocity products were computed. 
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2.2 J h i r d-Qrfls r ._C., ),<?.§ yrg.M.gM 

In view of the discussion given In the last part of the introduction 
section, the triple velocity fluctuation products, ?h'rjld be 

evaluated appropriately In a separating and reattaching flow region. This may 

“3 p 

be attained by formulating transport equations that solve for u , u v, 

5 ~3 

uv and v . The brief procedure of formulation and Its closure model Is 
given as follows. 

The general transport equation can be derived by suitable multiplication 
of the transport equation for a single fluctuating velocity u by unaveraged 
Reynolds stresses. Then, after time averaging of the resultant equation, the 
following is obtained. 


DUi Du, Du. Du.u.u. 
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This will lead us to the transport equation of u^u^ as 
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In closing Eq. (20), we need to make some assumptions. Because the first and 
second terms (I and II) on the RHS consist of mean strains, stresses, and 
triple products, they can be evaluated directly without making any correlation. 
The third term contains quadruple velocity fluctuation products, and these 
must be correlated In terms of known quantities such as the Reynolds stresses 
or mean strain rates. Mill lonshtchi kov 1 1 assumed that the quadruple 
velocity correlations are Gaussian and thus approximated them by using the 
product of the second-order correlations obtained by the formulae which are 
strictly fulfilled for the normal law. Then the quadruple velocity product 
can be decoupled as 

u i“jVi + + < 21 > 

and, hence, the sum of terms II and III is finally rearranged as 


II + III = 


Cu . u_ 
i si 


au,u k _____ 

+ v* 


au k u i 
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+ u, u — ! ~ J *1 
k i ax, J 


( 22 ) 


The last term (IV) can be correlated with the triple correlation as 

IV - - c s k VA (23) 

where the coefficient C $ needs to be determined by a computer optimization 
technique. Through our preliminary examination, it is found that the 
computations with C g values of between 1.0 and 10.0 result In acceptable 
levels of "u 1 u j u k . 


- 9 - 


With the approximation mentioned above, Eq. (20) Is finally written as 
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2.3 Other Models of Third-Order Closure 

The transport model of the third-order velocity products described in the 
preceding section was developed along with other existing models for the flow 
field of reattaching shear layer. In this section these other models are 
discussed. 

5 

Daly and Harlow obtained an algebraic expression with a simplification 
of such a transport equation of u^u^u^ as follows 


Vj u k - 
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1 2 

Hanjalic and Launder neglected the first term on the right-hand side 

of Eq. (24) based on their experimental measurements of axisymmetric flow in a 
1 3 

plane channel. Furthermore, upon neglecting the convective transport of 
Eq. (24), the following algebraic equation was obtained. 


au j\ 


u.u.u k = - 0.08 C 3 7 [u.u^ + U JU| 


!Vi 

ax. 


* Vi 


ax. 


(26) 


- 10 - 


Shir 


14 


proposed the following expression 


Vj u k ' - °- 04 C 3 « 
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Cormack et al. obtained algebraic expressions by approximating the 
experimentally determined profiles for u^Uj , c and with a 

polynomial with coefficients chosen to give a least-square fit to the data and 
discretized the cross-stream coordinate In each flow Into p grid intervals by 
using the most general model for the triple velocity correlation tensor as 
generated using the asymptotic approach of Lumley and Khajeh Nouri, 7 The 
model they obtained Is as follows: 


4|/* 
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(28) 


where 


a ij " Vj " 3 k4 1j 


(29) 


Out of the twenty parameters, a^, that they had started out with, 
they were able to determine and optimize the most significant ones as being 
a.|, <* 5 , ay and a 12 . The values of these parameters have been 
recommended for various kinds of flows along with the universal value 
applicable to most of the flows. 
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<*| - « 8.14 X fO * 
a K * - 1.72 X TO" 2 

( 30 ) 

Oy ss -4.80xl0 2 
°12 * ” 

In the original modes of Eqs. (25)-(28), the coefficient C 3 was chosen 
to be 1, but the best values for C 3 were Investigated for the reattaching 
shear flows by comparison with experimental data. The recommended values are 
given In Sec. 4.3. 


3. NUMERICAL MODEL 


3.1 Numerical Method 

The solution method of the transport equations described in the preceding 

1 6 

section Is the same finite volume method as that employed in TEACH Code, 

while the differencing scheme is the modified hybrid scheme of Amano^ 7 In 

which the combined mode of convection and diffusion is derived by expanding 

the analytical one-dimensional solution up to the fourth-order term. The cell 

structure for mean-velocity components is the staggered system in which the 

locations of the mean velocities U and V are a half-cell shifted in x- and y- 

~~2 

directions, respectively. All the normal Reynolds stresses (u^) are 
evaluated at the scalar node point along with P, k and c. However, the 
shear Reynolds stress (u7) is located at the southwest corner of the scalar 
cell. This is because the main driving strains for the shear stress are 
aU/ay and aV/ax which can easily be evaluated without any 
Interpolations (Fig. 1). 
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3.2 Boundary Conditions 


There are three different types of boundary conditions to be specified 
for the computation of a flow In the channel as shown In Fig. 2: Inlet, 
outlet, and wall boundary conditions. 

At the inlet all the quantities are specified according to the fully 
developed condition. At the outlet a contlnuatlve flow condition Is applied 
where gradients of flow properties In the flow direction are zero (Neumann 
conditions), l.e,, 3<(>/ax * 0, where 4> « U^, k, e, u^juj, etc. This outlet 
is located about 80-1 20H downstream from the step so that its Influence on the 
main flow region Is negligibly small. 

At the wall boundaries the velocities and turbulence quantities must be 
specified functionally according to the drag law or the law of the wall. For 
example, the tangential velocity can be expressed In terms of wall shear 
stress as a functional expression of the boundary condition coupled with the 
no-slip condition as 

F = tAA (31) 

where 6A = wall area of the cell. The velocity component normal to the wall 
is simply set as zero. The wall boundary values for k and c are determined 
by means of wall functions based on the assumption of a logarithmic near-wall 
velocity distribution which allows the wall-shear stress to be extracted from 
the "log law" and the value of velocity parallel to the wall to be computed 
along the grid line closest to the wall. Near-well effects on the turbulence 
structure, associated with the steep velocity variations, are also taken Into 
account by introducing appropriate modifications to the generation and dissi- 
pation of the turbulent energy and the energy dissipation rate for the finite 

1 7 

volume adjacent to the wall. 
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The boundary values for the Reynolds stresses are determined as 


u 2 » 1.21k 

v 2 - 0.24k (32) 

^uv » -0.24k + * ^ 

in the wall adjacent numerical cells. These boundary conditions are developed 
noting the following points: i) the coefficients in the normal stresses 
represent a consensus of several reported wall flows, and ii) the boundary 
condition for the shear stress Is given by the mean momentum equation. The 
details of the derivation of Eq. (32) are given in the Appendix. 

3.3 Computing Details 

Exploratory tests were made for different mesh sizes to investigate an 
optimum grid-independent state. Figure 3 represents a length of reattachment 
for different grid expansion factors in x-direction ranging from 1.01 to 1.05 
and for different numbers of grid points in x- and y- directions. The larger 
value of the expansion factor e x creates a finer mesh in the recirculating 
region while it results in a coarser mesh size in the downstream region far 
from the step. As shown in the figure, the variation In the reattachment 
length does not exceed 2 % between the medium value of 1.03 and the relatively 
fine value of 1.05. Also, the grid system almost reaches a grid-independent 
state around 40 grid points for both x- and y- directions. For most of the 
computations a 42 x 42 grid system is used with an expansion factor of 1.03 
in x-direction. However, it should be emphasized that the above mentioned 
test is made only for a prediction of a reattachment length. At this stage we 
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simply carried out most of the computations based on this test leaving the 
refinement of the grids for the Reynolds stresses or other turbulence 
quantities to the next progress report. 

In the computation of the Reynolds-stress equations (7), Initial stress 
values have to be specified properly so as not to cau>e numerical divergence 
since the system of equations Is relatively unstable. This Instability Is 
mainly due to an explicit form of the Reynolds stress equations which contains 
Reynolds stresses In most of the terms on the right-hand side. Thus a new 
value of iTjuJ Is comprted by Inserting previously computed ujuj Iteratively 
until it converges. This Instability problem was overcome by employing a 
three-pass procedure, that is, the computation was initially started with the 
standard k ~ c model. The Boussinesq viscosity model of the Reynolds 
stresses was Invoked when the maximum residual source of the mass and momentum 
became 5 percent of the total Inlet mas./ or momentum which was attained In 
about 50 iterations. The transport equations of the Reynolds stresses were 
brought In when this residual source became 3 percent which took another 20 
iterations. Complete convergence was assumed to have been attained when the 
relative maximum residual source of the Reynolds stresses was less than 1.7 x 
10" 3 which required typically another 320 Iterations with a total CPU time 
of about 50 minutes on a UNIVAC 1100 computer. 

4. RESULTS AND DISCUSSION 


4.1 Mean Velocity Profiles 

The flow in the channel shown in Fig. 2 was computed by using the models 
described In the preceding section, and the results were compared with 
available experimental data. 
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Figure 4 shows the velocity profiles at several different locations along 

the flow-stream for the channel step ratio Y Q /H « 3 and for Re e - 5000. 

Both the computations by the k ~ e model and the RSM were compared with 

1 8 

the experimental data of Seegmlller and Driver , at three streamwlse posi- 
tions (x/H * 2, 5.56, and 8.63), where the reattachment point was about x r /H 
* 5.0. Thus, x/H ■ 2 corresponds to the position In the recirculating 
region, x/H * 5.56 to the position near the reattachment point, and x/H * 8.63 
Is the position In the flow redeveloping region. It Is observed that 
agreement between the experimental data and computational results In the 
recirculating and reattachment regions Is not as good as that In the 
redeveloping region. This Is because the numerically predicted reattachment 
length, x r /H - 5, Is somewhat different from that determined by the 
experiment, x r /H = 7. Moreover, although the difference between the 
computations by the k ~ c model and the RSM is not remarkable, the RSM 
consistently gives slightly better results. 

Figure 5 shows the velocity profiles for the step ratio Y Q /H = 2 at 
three locations: x/H = 2.67, 5.33 and 9.78. The location x/H = 5.33 is very 
near the reattachment point since the computed x r /H is about 5.5. In 

addition, the computations are compared with the experimental data of Kim et 

* 

1 9 

al. where the trend of the prediction is quite similar to that in Fig. 4. 
However, in this case agreement between experiment and computation in the 
recirculating region is better than that found near the reattachment point. 

The disagreement near the reattachment point Is because there Is a discrepancy 
in the prediction of the reattachment length between experiment <x r /H ■ 7) 
and computation (x f /H = 5.5). 


- 16 - 


In Figs, 4 and 5 the computations by an algebraic stress model (ASM) made 
20 

by Slndlr are also compared. The best agreement for the recirculating 
region Is predicted by the ASM and for the redeveloping region by the RSM. 

As discussed above, the turbulence model Is dependent on the location of 
the flow field. In the recirculating region the RSM provides slightly better 
predictions near the wall. However, In the shear flow part, the ASM gives the 
best agreement with measurements. In the reattachment region the RSM shows 
the slowest velocity recovery performance while at the same time In the 
redeveloping region It shows the best prediction. Since the behavior of the 
mean velocity is closely related to the performance of Reynolds stresses, we 
will next examine the action of normal and shear stresses In the flow region 
across the streamwise flow. 

4 . ?, Reynolds Stresses 

Figures 6 and 7 show, respectively, the normal and shear stress profiles 

for the channel expansion of Y Q /H = 3, which correspond to the velocity 

1 8 

profiles In Fig. 4. The experimental data of Seegmiller and Driver are 
also shown in these figures. The normal stresses Increase rapidly with the 
distance from the bottom wall and reach their maxima at about y/H =1.0 in the 
recirculating region and at y/H = 0.5 - 0.7 in the reattachment and redevel- 
oping regions. The stresses then decrease to very small values at approxi- 
mately y/H = 1.3 - 1.5. It Is of interest to note that the y location of the 
peak moves upwards in the order of the k ~ c, RSM, and ASM. The level of 
normal stresses also increases In the same order as the peak location. A 
similar trend is also observed for the shear stress distributions except that 
they have negative values. 
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With regard to the comparison between measurement and computations, 
the predictions by the ASH seem to be better for the normal stresses while 
those by the RSM are closer to the experimental data for the shear stresses. 
The predicted levels of the Reynolds stresses by the k ~ c model are 
always too low. 

Observation of these Reynolds stresses Indicates that the performance of 
the Reynolds-stress prediction does not necessarily accord with the compu- 
tation of the mean velocity profiles which are shown In Fig. 4. At x/H « 2, 
the RSM gives better agreement, but the ASM Is superior for normal stresses in 
the redeveloping region. Generally, we see the success of the Reynolds-stress 
closure In the prediction of the shear-stress distributions. 

4.3 Triple Velocity Products 

Computations for the triple velocity products were performed for two sets 

of flow conditions corresponding to the experimental data of Driver and 

21 22 
Seegmiller and Chandrsuda and Bradshaw. 

Figures 8, 9 and 10 show the profiles for the normalized triple velocity 

products (uuv + vvv)/Uj N as a function of nondimensional distance from 

the step y/H at three different locations of x/H = 4.0, 6.0 and 8.0. The 

computation was performed using a 42 x 42 mesh with cell expansion factors of 

1.03 and 0.98 in the axial and transverse directions, respectively. Triple 

velocity products were obtained by the four correlations proposed by Daly and 

5 1214 15 

Harlow , Hanjalic and Launder , Shir , and Cormack et al. , and are 

21 

compared with the experimental data of Driver and Seegmiller for a 
backward-facing step with a step ratio, Y q /H = 8. 

The models given in Eqs. (25)~(28) were used to compute the triple 
velocity products after all the converged solutions of the Reynolds stresses 
were obtained. 
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All these models for the triple velocity correlations were Initially 
developed to satisfy the conditions of relatively simple types of flow such as 
free jets, wakes, boundary layers, etc. Therefore It Is necessary to modify 
these models for the present reattaching and recirculating flows. The 
coefficient, C g , In Eqs. (25>~( 28) was examined In order to determine 
optimum results of triple velocity products throughout the solution domain. 

Upon comparing the coefficients C 3 In these correlations with the 
experimental data, It was found that the value of C g differed by a factor of 
1 to 50 among different correlations. The determination of this empirical 
multiplication factor, C 3 , Is based upon a simple averaging method wherein 
the peak maxima and minima are measured and compared with the experimental 
data and factors were thus obtained. 

The factors developed and used are given In Table 1. 

Table 1. Factors used In different models 


Model 

c 3 

Shir 

53.13 

Hanjalic and Launder 

6.27 

Cormack et al . 

5.67 

Oaly and Harlow 

1.06 


The general trend prescribed by these correlations is very similar to the 

experimental data as seen in the Figs. 8, 9 and 10. 

Figures 11, 12 and 13 compare the experimental data of Chandrsuda and 
22 

Bradshaw for a backward-facing step with a step ratio Y Q /H * 2.5 with 
the triple velocity products obtained by the aforementioned four correlations 
at three different locations of x/H = 5.4, 8.4, and 10.3, respectively. These 
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computations were performed using a 52 x 52 mesh with cell expansion and 
contraction factors of 1.015 and 0,990 In the axial and transverse directions, 
respectively. The computational technique was similar to the one used to 
compute the data for Y Q /H 8 except for the fact that the convergence was 
assumed to have been completed when the maximum residual source was 1.5 
percent (which is Independent of the source term for Reynolds stresses), which 
took about 620 Iterations total. In this computation about 20U Iterations 
were performed for solution of the transport equations of the Reynolds 
stresses. Total CPU time on a UNIVAC 1100 computer was about 90 minutes. 

Figure 11 compares the variation of vvv/uJ N as a function of y/H at 
the location x/H * 5.4 with different correlations and the experimental data. 
Figure 12 compares tie variation of »uv/uJ N as a function of y/H at the 

- o 

location x/H « 8.4, and Fig. 13 compares the variation of uuv/Uj N at the 
location x/H * 10.3. There Is generally good agreement between the results of 
the correlations and the experimental data as seen in these figures. 

It Is rather premature to discuss the relative merits and demerits of 
each correlation of triple velocity products at this stage. However, It has 
been seen that most of the correlations presented here may be Improved with 
only minor modifications. 


5. SUMMARIZING REMARKS 

Through the first six month period, a hybrid model of the Reynolds-stress 
closure was developed. This model was tested for various sizes of step flow, 
and the computed Reynolds-stress behavior was compared with experimental data 
as well as with other simpler turbulence models such as the algebraic stress 
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model and the standard k ~ c model. Although the results are not 
drastically Improved with the RSM, It Is shown that overall agreement with 
experimentally measured values for both the mean velocity and the Reynolds 
stresses Is better than other models. Moreover, the application of the RSM to 
flows In geometrically more complicated shapes Is promising. It Is also noted 
that the k ~ c model or the ASM cannot be used for swh complex flow 
predictions due to the fact that the k ~ e model Is limited only to a 
quasl-lsotroplc flow and the ASM Is applicable only when the convective effect 
of the Reynolds stress is not Important. 

Secondly, the third-order closure model was reviewed and transport 
equations for the triple velocity correlation were developed and Implemented 
In a numerical code to evaluate the behavior of the triple velocity products 
in various regions of the flow field including recirculating, reattaching t ind 
redeveloping How domains. These four existing models of algebraic equations 
for the trlpte velocity products were also employed In a program and a 
modification of the empirical coefficients was made to improve the prediction 
for the reattaching shear flows. 

Probably the most significant fact the Investigators have learned through 
this project Is that both convection and diffusion effects of the triple 
velocity products become very important in the reattaching shear layer which 
has been Ignored by most of the researchers. It Is our aim to pursue the 
effects of these rates In the third-order closure model. In the next stage of 
this project, the transport model of the third-order closure will be 
Implemented and tested for universal usage. 


6, WORK IN PROGRESS 


6j QMigK m m 

Although all the computations made In this report are done by using the 
modified hybrid scheme specified in Sec. 3.1, the authors are seeking a better 
discretization method for computing the transport equations used In this 
project. 

When convection of a transport equation does not predominate over 
diffusion, the exclusive use of central differencing in space undoubtedly 
offers a satisfactory compromise between the requirements of accuracy and 
computational economy. However, when convection Is dominant, central 
differencing for convection results in Instability. To avoid this 
Instability, a false diffusion Is frequently brought In at an expense of 
accuracy. 

The reduction of false diffusion In the computation of steady-state 

recirculating flows at no significant loss of computational economy has been 

one of the objects of this project. Some of the schemes proposed by several 

researchers have been reviewed. The "skew-upwind-differencing scheme" of 

Raithby yields a certain amount of reduction in skewness error by 

partially simulating an upwind discretization in a stream-line coordinate 

system. By using Raithby's method, skewness errors are totaly absent although 

24 

it Is only first-order accurate. The method of Leonard' Is of a more 
fundamental nature, designed to eliminate false diffusion altogether by using 
a convective manner, i.e., quadratic, upstream-weighted Interpolation for the 
convection terms (QUICK). 

The original work of the QUICK method demonstrated by Leonard was for 
one-dimensional and for uniform-grid system only. Thus this method needs to 
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be tailored for more general two-dimensional and non uniform grid systems. 
The procedure for this Is given as follows. 

Consider the problem of approximating * w In Fig. 14. Depending on 
the sign of U either points WW, W and P or points W, P and E are used to 

W 

fit a quadratic curve for Interpolating the value of * at w. A length 
i 3 Is now defined as the distance from the west face to either the 
WW-poInt or E-point depending on the direction of U . It can be written as 

WW 


6 



A *P£ + *2 U w <0 

. " (AX WW-W * V U W > 0 


(33) 


A quadratic fit then yields for 4 > w as 


4*, 


w 


r ■*« 




where 


4* 


<> WW U M > ° 

*E U w < ° 


and 


B 1 = V3 (< 2 ' 4 3) 
B ? = 4,43(4, + 4 j) 

B 3 = 4,4 ? (4, ♦ 4 ? ) 
B - B, - Bj + B3 


The above formulations are Implemented for a two-dimensional system. 


(34) 


(35) 


(36) 
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So far the Implementation of QUICK Is done only for the momentum 
equations— It Is still being tested for the rest of the equations. Therefore 
it seems to be premature to conclude any advantages and disadvantages 
connected with its use. However, we feel that there would be only small 
differences In the results of the third-order closure model no matter which 
system is employed. As Is noted in most of the transport equations, 
generations and dissipations are mainly evaluated through the source terms of 
the discretization equations not through the convection-diffusion 
discretization. 

6.2 Consolidation of the RSM 

The work described in the preceding sections is not a complete model of 
the full Reynolds-stress equations. The diffusion terms of the momentum 
(Eq. (2)) and the dissipation equations (Eq. (£)) still contain the Boussinesq 
viscosity approximation simply because the solutions become numerically 
unstable when we employ the Reynolds stresses in these diffusion terms. It 
was also mentioned that the third-order closure was employed only after all 
the converged solutions for the Reynolds stresses had been obtained. The 
following briefly summarizes the work currently being conducted to accomplish 
the ceding of the full Reynolds-stress model. 

i) Diffusion Rate of Momentum Equations. 

The form of Eq. (2) was converted into the following form: 


3P 


.j <' u iV - - S7 + 


aU i 


au. 


■j [V S + 


ax* ax. ■ ax. LH 'ax. ’ ax/ au./ax* ax 


k i 


pu.u. au, 

A _ I_J 1/1 

' All /av Av J 


{ 37 ) 




Equation (37) can be cast into the standard form as 
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(38) 


Sj <' u iV 


4£- + -2- r (u _ -fVi-, 

ax, axj Uw au,/axj' 


au 1 a au 

3X j^ 4 3Xj 


) 


In solving Eq. (38) numerically, the first and the third terms on the 
right-hand side are evaluated In the source term of the discretization 
equation, whereas the second term Is treated as a diffusion term where the 
effective viscosity becomes 


eff 


pu^i 

at^/axj 


(39) 


This treatment Is found to cause less Instability and the convergence process 
Is Improved. 

2) Diffusion Rate of c-Eouatlon . 

The dissipation equation (Eq. (5)) which Is used to solve the RSM, Is a 
second-order model except for the diffusion term which is again evaluated by 
using the Boussinesq viscosity model. The complete second-order c-equation 
can be obtained by replacing the last term in Eq. (5) by the following 


0(c) = 


3x j 


C(m + C 


ek 

c 


. dc 

u i u j^ax. 


(40) 


The mechanism and significance of both primary and secondary diffusion terms 
are now being investigated. 

3) Merging the Triple Velocity Products in the RSM. 

As mentioned earlier in this report, the behavior of the third-order 
tensor has a strong effect particularly near the reattaching flow region. 

This will considerably influence the diffusive action of the RSM. This is the 
main object of the modeling of the third-order closure. 
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For the same reason as stated In case (1) of this section, the form of 
Eq. (11) Is recommended to be modified to a formal diffusion term such as 
those given by Eq. (18). Consequently, Eq. (11) was transformed Into the 
following 



u 1 u 1 u k 

3 <Vj )/ax m 


au^ 


ax 


m 


( 41 ) 


The computation of Eq. (41) Is now being undertaken In parallel with the 
development of the third-order closure model. 
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7. APPENDIX 


7.1 Near-Wall Values for Reynolds Stresses 

At the wall the boundary conditions for the Reynolds stresses are 
determined as follows. First, the boundary condition on uv Is given by the 
mean momentum equation (Eq. (2)), Near the wall where convection of the 
momentum is negligibly small, Eq. (2) yields the following form. 

0 * " d? + dy (T w ' <*•’> 

It Is also known from the Law of the Wall that the wall shear stress may be 
given as 

U* - (A. 2) 

Thus, by integrating Eq. (A.l) and after substituting Eq. (A. 2), we obtain the 
following. 


- uv = - 



T 


y /d£ 

^dx 


(A. 3) 


The determination of the normal components of the Reynolds stresses is 
based on an observation of experimental data. A consensus of several of the 
most -thoroughly documented wall flows [25-28] is represented as 

u 2 = 4.5 U 2 

T 

V 2 = 1.0 U 2 (A. 4) 

~2 2 
W* = 2.0 IT 
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By using the correlation of the turbulence energy, 
velocity, U^, In the wall proximity region, 

k - 3.7 U 2 

the va^es of Eqs. (A. 3) and (A. 4) become 

- uv - -0.24k + ^ 

p dx 

- u 2 = 1.21k 

- v 2 - 0.24k 

- w 2 = 0.55k 

It should be noted that the boundary condition given by 
the definition of k, l.e., 

k = (u 2 + v 2 + w 2 )/2 


, and the friction 


(A. 5) 


(A. 6) 


. (A. 6) satisfies 


(A. 7) 


- 28 - 


8 . 


REFERENCES 


1. Bradshaw, P. and Wong, F. Y. F., "The Reattachment and Relaxation of a 
Turbulent Shear Layer," Journal of Fluid Mechanics . Vol. 52, Part 1, 

1972, pp. 113-135. 

2. Smyth, R., "Turbulent Flow Over a Plane Symmetric Sudden Expansion," 
Journal of Fluids Engineering . Vol. 101, No. 3, 1979, pp. 348-353. 

3. Eaton, J. K. and Johnston, J. P., "An Evaluation of Data for the 
Backward-Facing Step Flow," Stanford Conferences on Complex Turbulent 
Flows, 1980. 

4. Rotta, J. C., "Statlstlsche Theorle Nlchthomogener Turbulenz," Z. Phvs. . 
Vol. 129, 1951, p. 547. 

5. Daly, 8. J. and Harlow, F. H., "Transport Equations of Turbulence," 
Physics of Fluids . Vol. 13, No. 11, 1970, pp. 2634-2649. 

6. Naot, D., Shavlt, A. and Wolfshtein, M., "Two-point Correlation Model and 
the Redistribution of Reynolds Stress," Physics of Fluids . Vol. 16, 

No. 6, June 1973, pp. 738-743, 

7. Lumley, J. L. and Khajeh Nourl, B., "Computational Modeling of Turbulent 
Transport," A dvances In Geophysics . Vol. 18A, 1974, pp. 169-192. 

8. Launder, B. E., Reece, 6. J, and Rodl, W., "Progress in the Development 
of a Reynolds-stress Turbulence Closure," Journal of Fluid Mechanics . 

Vol. 68, 1975, pp. 537-566. 

9. Rodl, W., "The Prediction of Free Boundary Layers by Use of a Two- 
Equation Model of Turbulence," Ph.D. Thesis, University of London, Dec. 
1972. 

10. Amano, R. S. and Goel, P., "A Numerical Study of a Separating and Reat- 
taching Flow by Using Reynolds-stress Turbulence Closure," Numerical Heat 
Transfer . Vol. 7, No. 3, 1984, pp. 343-357. 

11. Millionshtchikov, M. D., "On the Theory of Homogeneous Isotropic 
Turbulence," C. R. Acad. Sci. S.S.S.R. . Vol. 32, 1941, pp. 615-619. 

12. Hanjalic, K. and Launder, B. E., "A Reynolds Stress Model of Turbulence 
and Its Application to Thin Shear Flows," Journal of Fluid Mechanics . 

Vol. 52, Part 4, 1972, pp. 609-638. 

13. Hanjalic, K. and Launder, B. E., "Fully Developed Asymmetric Flow in a 
Plane Channel," Journal of Fluid Mechanics . Vol. 51, Part 2, 1972, 

pp. 301-335. 

14. Shir, C. C., "A Preliminary Numerical Study of Atmospheric Turbulent 
Flows in the Idealized Planetary Boundary Layer," Journal of Atmospheric 
Science . Vol. 30, 1973, pp. 1327. 


29 - 


15. Cormack, 0. E. , Leal, L. G. and Seinfeld, J. H., "An Evaluation of Mean 

Reynolds Stress Turbulence Models: The Triple Velocity Correlation," 

Transactions of the ASME . Journal of Fluids Engineering . Vol. 100, 

March 1978, pp. 47-54. 

16. Gosman, A. 0. and Ideriah, F. J. K., "TEACH-T: A General Computer Pro- 
gram for Two-Dimensional, Turbulent, Recirculating Flows," Mechanical 
Engineering Department Report, June 1976, Imperial College, London. 

17. Amano, R. S., "Development of Turbulent Near-Wall Model and Its Applica- 
tion to Separated and Reattached Flows," Numerical Heat Transfer . Vol. 7, 
No. 1, 1984, pp. 59-75. 

18. Seegmlller, H. L. and Driver, D. M., Private Communications, 1984. 

19. Kim, J., Kline, S. J. and Johnston, J. P., "Investigation of Separation 

and Reattachment of a Turbulent Shear Layer: Flow Over a Backward Facing 

Step," Report MD-37, Thermosciences Division* Department of Mechanical 
Engineering, Stanford University, 1978. 

20. Sindlr, M. M., "Effects of Expansion Ratio on the Calculation of 

Parallel-Walled Backward-Facing Step Flows: Comparison of Four Models of 

Turbulence," ASME Paper 83-FE-10, 1983. 

21. Driver, D. M. and Seegmlller, H. L., "Features of a Reattaching Turbulent 
Shear Layer Subject to an Adverse Pressure Gradient," AIAA Paper 82-1029, 
1982. 

22. Chandrsuda, C. and Bradshaw, P., "Turbulence Structure of a Reattaching 
Mixing Layer," Journal of Fluid Mechanics . Vol. 110, 1981, pp. 171-194. 

23. Raithby, G. D., "Skew Upwind Differencing Schemes for Problems Involving 
Fluid Flow," Computer Methods in Applied Mechanics and Engineering . 

Vol. 9, 1976, p. 153. 

24. Leonard, B. P., "A Stable and Accurate Convective Modelling Procedure 
Based on Quadratic Upstream Interpolation," Computer Methods in Applied 
Mechanics and Engineering . Vol. 19, 1979, pp. 59-98. 

25. Laufer, J., "Investigation of Turbulent Flow in a Two-Dimensional 
Channel," NACA Report . No. 1053, 1951. 

26. Klebanoff, P. S., "Characteristics of Turbulence in a Boundary Layer with 
Zero Pressure Gradient," NACA Report . No. 1247, 1955. 

27. Chevray, R. and Kovasznay, L. S. G., "Turbulence Measurements in the Wake 
of a Thin Flat Plate." AIAA Journal . Vol. 7, 1969, p. 1641. 

28. Champagne, F. H., Harris, V. G. and Corrsin, S., "Experiments on Nearly 
Homogeneous Shear Flow," Journal of Fluid Mechanics . Vol. 41, 1970, p, 81. 


- 30 - 


9. FIGURES 


Figure 1 
Figure 2 
Figure 3 
Figure 4 
Figure 5 
Figure 6 
Figure 7 
Figure 8 

Figure 9 

Figure 10 

Figure 11 

Figure 12 

Figure 13 


Cell structure. 

Flow domain. 

Reattachment length for different grid size. 

Mean streamwlse velocity profile. 

Mean streamwlse velocity profile. 

Reynolds normal stress profile. 

Reynolds shear stress profile. 

Triple velocity products at x/H « 4. 

(comparison with experiment of Driver and Seegmlller) 

Triple velocity products at x/H ■ 6. 

(comparison with experiment of Driver and Seegmlller) 

Triple velocity products at x/H ** 8. 

(comparison with experiment of Driver and Seegmlller) 

Triple velocity products at x/H * 5.4. 

(comparison with experiment of Chandrsuda and Bradshaw) 

Triple velocity products at x/H = 8.4. 

(comparison with experiment of Chandrsuda and Bradshaw) 

Triple velocity products at x/H = 10.3. 

(comparison with experiment of Chandrsuda and Bradshaw) 


Figure 14 


QUICK for non-uniform grid. 
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Figure 3 Reattachment length for different grid size. 
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Moan streanwise velocity profile. 
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Figure 8 Triple velocity products at x/H = 4. 

(comparison with experiment of Driver and Seegmiller) 




Figure 9 Triple velocity products at x/H = 6. 

(comparison with experiment cf Oriver and Seegmiller) 
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Figure 10 Triple velocity products at x/H = 8. 

(comparison with experiment of Driver and Seegmiller) 




Figure 11 Triple velocity products at x/H = 5.4. 

(comparison with experiment of Chandrsuda and Bradshaw) 
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Figure 12 Triple velocity products at x/H - 8.4. 

(comparison with experiment of Chandrsuda and Bradshaw) 
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Figure 13 Triple velocity products at x/H = 10.3. 

(comparison with experiment of Chandrsuda and Bradshaw) 





